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ABSTRACT 

The high quality light curves of Kepler space telescope make it possible to analyze the optical 
variability of AGNs with an unprecedented time resolution. Studying the asymmetry in variations 
could give independent constraints on the physical models for AGN variability. In this paper, we 
use Kepler observations of 19 sources to perform analyses on the variability asymmetry of AGNs. We 
apply smoothing-correction to light curves to deduct the bias to high frequency variability asymmetry, 
caused by long term variations which are poorly sampled due to the limited length of light curves. 

A parameter /3 based on structure functions is introduced to quantitively describe the asymmetry 
and its uncertainty is measured using extensive Monte-Carlo simulations. Individual sources show no 
evidence of asymmetry at timescales of 1 ~ 20 days and there is not a general trend toward positive 
or negative asymmetry over the whole sample. Stacking data of all 19 AGNs, we derive averaged /3 
of 0.00±0.03 and -0.02±0.04 over timescales of 1^5 days and 5~20 days, respectively, statistically 
consistent with zero. Quasars and Seyfert galaxies show similar asymmetry parameters. Our results 
indicate that short term optical variations in AGNs are highly symmetric. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: Seyfert — quasars: general 


1. INTRODUCTION 

Aperiodic optical/ultraviolet variability is a significant 
property of Active Galactic Nuclei (AGNs), but its phys¬ 
ical origin is still unclear. Variability asymmetry de¬ 
scribes whether the light curve favors a shape of rapid 
rise and gradual decay, i.e., positive asymmetry, or a 
shape of gradual rise and rapid decay, i.e., negative asym¬ 
metry. There are only a few observational/theoretical 
works in literat ure studying the variab ility asymme¬ 
try of AGNs. iKawaguchi et al.l (|1998ll introduced a 
structure function approach to estimate the variability 
asymmetry of AGN light curves. They adopted two 
structure functions, i.e., SFic^r) and SF^dT), which 
only include pair epochs with increasing and decreas¬ 
ing flux, respectively. Possible physical models which 
could produce asymmetry in variations were also dis¬ 
cussed iii__this_work. Throu gh Monte Garlo simula¬ 
tions, IKawaguchi et al.l (II998I 1 showed that the disk in¬ 
stability model produces SFic{T) < SFdc{T) or negative 
asymmetry, while the starburst model, which attributes 
the optical variations to random superposition of super¬ 
novae in the nuclear starburst region, yields a contrary 
asymmetry, i.e ., SFic{T) > SFddT). Shortly thereafter, 
iHawkin^ (|2002f) added that the micro-lensing model pre¬ 
dicts SFicir) = S FdrjT) or no asy mmetry statistically. 

Observationally, iHawkinj (I2002f) tested the aforemen¬ 
tioned three models using long term optical light curves 
of 401 quasars and 45 Seyfert galaxies. They found 
that Seyfert galaxy NGG 5548 appears negative asym¬ 
metric variations on timescales of 10 ~ 70 days, while 
for quasars the variat ions are symmetric o n timescales 
of a year and longer. Ide Vries et al.l (|2005ll studied the 
long term variations of a large sample of 4I39I quasars 
with SDSS and historic photometry, and showed that 
on timescales of years the qua s ar var iations behave pos¬ 
itive asymmetry. I Bauer et al.l (1200911 analyzed the opti¬ 
cal variability of nearly 23000 quasars in the Palomar- 


QUEST Survey and found no evidence of any asym- 
metry in variabi lity over ^ 10 days to several years. 
IVoevodkiiil (|2010ll computed the structure functions for 
7562 quasars from SDSS Stripe 82 and detected a signif¬ 
icant negative asymmetry on timescales longer than 300 
days. All the analyses above were based on st r uctur e 
function me thod introduced by Kawaguchi et al.l (1199811 . 
Meanwhile, IGiveon et al.l (|1999 1 made use of a different 
method of calculating the difference between medians of 
brightening phases and fading phases in the light curves 
of 42 PG quasars and reported a negative asymmetry in 
variations. To summarize, negative asymmetry is favored 
by more works, but inconsistencies exist among the few 
observational studies in literature. 

The Kepler s pace telescope was designed to search 
for exoplanets (IBorucki et al.l I2010H and it could pro¬ 
duce nearly continuous optical light curves for the targets 
within its field of view, including AGNs. There are a cou¬ 
ple of works in literature reporting t he high frequency 
variab ility analyses of Kepler AGNs. IMushotzkv et all 
(|2011ll calculated the power spectral density (PSD) func¬ 
tions for four Seyfert 1 galaxies in the Kepler field 
and obtained best-fit PSD power-law slopes of —2.6 ~ 
—3.3, considerably steeper than those of quasars at 
timescales of months to years which could be described 
as d amped random walk process (e.g., IKellv et Ml 
200911 . ComDl e mente d with ground based observations, 
Garini &: Rvlel (I20I2I1 presented a furthe r analysis of 
the Kepler light curve of Zw 229-15. IWehrle et al.l 
(I20I3I1 and iRevalski et al.l (|2013ll reported four radio- 
loud Kepler AGNs and calculated the PSD functions 
using light cu r ves st itched by a normalization method. 
lEdelson et al.l (|2013ll reported the variability analyses of 
a BL Lac object, i.e., W2R 1926-1-42, in the Kepler field. 

Meanwhile, the Kepler data provide the first ever op¬ 
portunity to study the short term variability asymmetry 
of AGNs, which is the aim of this work. We describe the 
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adopted Kepler AGN sample and light curve stitching 
process in Section [2j The methods of asymmetry anal¬ 
ysis and smoothing-correction are introduced in Section 
m In Section |4] we report our major results. Discussion 
is given in Section [5] and conclusions in Section [6l 

2. DATA REDUCTION 

2.1. Kepler AGN Source 

There are only a few cataloged AGNs in the Ke¬ 
pler’s field of view. The 19 AGNs we used in this 
work an d listed in Tabl e H] a re collected from lit¬ 
erature (iMushotzkv et all 12^1 H lEdelson et al.l 120121 : 
IWehrle et al.ll201.1l l. Following lHawkin^ (j2002f l. we divide 
the AGNs into two subsamples, including 9 quasars with 
Mj < —23.6 and 9 Seyfert galaxies with Mj > —23.6. 
The rest source (W2R 1926-1-42) in the table is a BL Lac 
object. 

We adopt the SAP FLUX in the Kepler light curve 
files, since PDCSAP FLUX (calibrated for systematic ef¬ 
fects, e.g., pointing and focus cha nges) may have mask ed 
out the intrinsic AGN variations (jCarini fc R,vlel[2012[l . 

2.2. Stitching Light Curves 

As a result of the Differential Velocity Aberration effect 
of Kepler telescope, the target flux is continuously redis¬ 
tributed among neighboring pixels, appearing as artificial 
long term variations in electron counts within the fixed 
optimal aperture a nd discontinuous light curves between 
adjacent quarters llKeol er Instr ument HandbookI 120091 
iKinemuchi et al.l l2012f) . iCarini fc Rvlel ( 201211 rescaled 
the Kepler light curves with coordinated ground based 
observations to connect di fferent quarters, howeve r, it is 
not a general approach. IKinemuchi et al.l (j2012[l gave 
another method, which increases the number of pixels 
in target mask for photometry with the PyKE tasks 
kepmask and kepextract. With reset target mask, a 
new light curve can be extracted from the counterpart 
Target Pixel File. The later method could reduce tar¬ 
get flux losses out of the aperture, but it might also get 
contaminat ion from nearby sources . 

Following iKinemuchi et ahl l)2012ll . we perform the pix¬ 
els re-extraction process as illustrated in Figure [1] Since 
not all the sources have enough observational quarters 
and clear surrounding and, during the process, enough 
‘halo’ pixels surrounding the target are needed, only four 
sources, i.e., Zw 229-15, W2 1925-b50, W2R 1904-b37, 
and GGRaBS J1918- I-493 7, can be stitched. The results 
are shown in Section l4.ll 

3. ANALYSIS METHODS 

3.1. Structure Function and Asymmetry Parameter 

The general definition of structure fun c tion and its 
properties were given by iSimonetti et al.1 (jl985[l . The 
first-order structure function SF{t) is defined as 

E + ^) - (1) 

A(t) ^ 

where N(t) is the number of data pairs for a certain 
time lag r, and / the observed source flusQ- To quantify 

^ Computing the structure function and then the asymmetry 
parameter using magnitude does not alter the results presented in 
this work. 


Zw 229-15 
13 quarters 



cyclical mask 


Figure 1. The original and stitched light curves with bin size of 
30 minutes (hereafter the same) of Zw 229-15. Vertical axis denotes 
the flux (count rate). The original data is shown in red, while the 
re-extracted light curves are shown in yellow, green, blue, and pink. 
Four reset target masks are adopted (plotted with different colors 
in the light curves) and duplicately used for four quarters in each 
year. The four masks are shown at the bottom for demonstration. 


variability asy mmetry, an asymmetry parameter /3 (t) is 
defined as Isee IKawaguchi et al.lll998ll 


/3(t) 


SF^c{t) - SFdcir) 
SFtotir) 


( 2 ) 


where the suffix ic and dc denote data pairs with in¬ 
creasing and decreasing flux, respectively, and tot for the 
total data pairs. /3(t) quantifies the normalized differ¬ 
ence between the increasing and decreasing variability. 
A positive /3(t) represents that the light curve favors a 
shape of rapid rise and gradual decay, i.e., positive asym¬ 
metry, while a negative /3 (t) characterizes the opposite 
situation. 


3.2. Smoothing-Correcttion 

Since the light curves have limited duration, the long 
term variations are poorly sampled and will yield unphys¬ 
ical and large scatter in the measurement of the asym¬ 
metry parameter. In other words, for individual light 
curves, asymmetry analyses could only be performed to 
short term variations (the ratio of the duration of the 
light curve to the concerned timescale of variations ^ 
1), which have been well sampled. Furthermore, the 
long term trend in the light curves could also introduce 
biases into the short term asymmetry analyses. Such ef¬ 
fect can be reduced using smoothing-correction method, 
in which the corrected light curve is produced by the 
inverse Fourier transform of the power spectrum with 
the low frequency part set to zero. The cutoff frequency 
in the power spectrum is set to 2x10“^ Hz, which corre¬ 
sponds to about 60 days and is shorter than the duration 
of light curves for each quarter of all sources. An exam¬ 
ple is illustrated in Figure [5] and Section S] presents the 
results of asymmetry analysis for both the original and 
smoothing-corrected light curves. 

3.3. Uncertainty Estimation 
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Figure 2. An example shows the smoothing-correction applied 
to the tenth quarter of Zw 229-15. The left two panels show 
light curves with the vertical axis denoting the flux (count rate), 
where the original (grey) and smoothing-corrected high-frequency 
(red) light curves are shown in the upper and lower panels, respec¬ 
tively. The low-frequency curve (green) in the upper-left panel is 
the difference between the original and the high frequency curves. 
The right two panels show the corresponding structure functions. 
SFic{T) and SF^ci'^) ^re shown in orange and blue, separately, and 
the total SFtot{T) in grey. Strong positive asymmetry, i.e., SFic{r) 
> SFdciT), is visible at timescales above a couple of days in upper- 
right panel, due to the long term increasing trend seen in the origi¬ 
nal light curve. Weak positive asymmetry is also visible at shorter 
timescales in the original structure function (upper-right panel), 
but disappears after smoothing-correction (lower-right panel). 

Estimating the u ncertainties in SF(t) and /3(r) is not 
straightforward fsee lEmmanoulopoulos et al.ll201Cll) . The 
[f{t + t) — /(r)] series (in Equation [T]) of a single light 
curve are not mutually independent, thus the uncertain¬ 
ties in SF{t) and /3(t) at given r would be significantly 
underestimated with the standard deviation definition. 
Furthermore, the values of SF{t) and /3(r) at differ¬ 
ent times cales r are not independent eith er. Therefore, 
following lEmmanoulopoulos et al.l (j201fll l. we estimate 
the uncertainties through extensive Monte-Carlos sim¬ 
ulations. 

We adopt the algorithm of iTimmer fc Koenid (|1995f) to 
generate artificial light curves for given power spectrum. 
In order to make the simulated light curves possessing 
similar power spectrum shape and noise level with the 
observed ones, we use the original periodogram, calcu¬ 
lated from the original observed (and not the smoothing- 
corrected) light curves, as the input spectrum, instead 
of using power spectrum with a specific fixed power- 
law slope. The observed light curve is end-matched be¬ 
fore periodogram calculation to reduce contamination 
caused by t he mismatch between beginning and end 
points (e.g.. iMushotzkv et al.ll20liL iWehrle et al.ll20ll : 
lEdelson et al.l I2014D PI. To simulate the effect of red 
noise leak to the short light curves, the artificial light 
curves should be much longer than the observed ones, 
by extending the input power spectrum to lower fre- 

^ Without end-matching the high frequency components of 
the power spectrum will be spuriously enhan ced, then the 
spectrum slope will be flatter than reality IIFougerd 119851 : 
IGonzalez-Martm & VaughanI [20121 ). The simulated light curves, 
based on such input power spectrum, will be significantly biased 
against the observed one. Note that in this work the end-matching 
was only adopted to calculate the power spectrum. No end¬ 
matching is performed for calculations of structure function and 
asymmetry parameter. 
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TIME (BJD-2454833) 


Figure 3. The stitched light curves of Zw 229-15. The verti¬ 
cal axis denotes the flux (count rate). The original (grey) and 
smoothing-corrected high-frequency (red) light curves are shown 
in the upper and lower panels, respectively. The low-frequency 
curve (green) in the upper panel is the difference between the orig¬ 
inal and the high frequency curves. Blue dotted lines represent 
boundaries of intervals of different quarters. 

quencies (e.g., lUttlev ef al.l l2002t IVatierhan et al.l 120031 
lEmmanoulopoulos et al.ll2010H . We fit the low frequency 
part of the original power spectrum to measure the low 
frequency extending slope. If the best-fitted slope is 
steeper than -2, we adopt a fixed value of -2, as seen in 
the observed (lower frequency c omparing with K epler) 
power spectrum of quasars (e.g., iKellv et all 120091 ) 0 . 

For each observed Kepler light curve, we generate a 
single, 3000 times longer light curve, which is then split 
into 3000 segments. We then randomly select 1000 seg¬ 
ments of them to calculate the corresponding SF{t) and 
/3(t), following exactly the same procedures as we applied 
on the observed light curves, and take their scatters as 
the uncertainties of the observed SF{t) and /3(t), re¬ 
spectively. Note that during the simulations, it has been 
assumed that the variations are intrinsically symmetric, 
thus the output mean of simulated /3(t) equals zero, and 
its scatter represents the uncertainty of the observed /3(t) 
if not severely different from zero. 

4. VARIABILITY ASYMMETRY RESULTS 

4.1. Results of Quarterly-Stitched Light Curves 

In this section, we present the results from asymmetry 
analysis of the stitched light curves, constructed for the 
four sources listed in Section 12.21 (i.e., Zw 229-15, W2 
1925-b50, W2R 1904-b37, and CGRaBS J1918-b4937). 
The stitched light curve of Zw 229-15, spanning ^ 3.3 
years with a time resolution of 30 minutes, is presented 

® Note that adopting steeper lower frequency power spectrum 
slope would yield slightly larger scatter in the derived SF{t) and 
/3 (t), due to stronger red noise leak. However, the asymmetric 
scatter of smoothing-correction data, which we mostly care about, 
will not be influenced by the enhanced low frequency leak. 
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Figure 4. SF(t) of stitched light curves of Zw 229-15. SFic{r) 
and SFdci'^) are shown in orange and blue, separately, and the 
total SFtot{T) in grey; shadow areas represent 2cr error range of 
SFtotiT). The top panel is for the original data and the bottom 
for the smoothing-corrected data. 

in Figure [31 Zw 229-15 is a narrow-line Seyfert 1 galaxy, 
observed by Kepler during 13 quarters (Q4^Q16, the 
most among Kepler AGNs). 

Figure SI shows structure functions of both the original 
and smoothing-corrected light curves of Zw 229-15. The 
original SF{t) shows a rise towards longer timescales 
with a gradually decreasing slope. At timescales above 
200 days, SFidr) is significantly larger than SF^cir), 
consistent with the general increasing trend in the orig¬ 
inal light curve. At very short timescales, the SF{t) 
reaches a plateau due to the photometric uncertainty of 
the light curve, with value of v^cTm, where am is the 
photometric uncertainty. On the long term end, SF(t) 
converges to '/2ai, where ai is the standard deviation 
of whole light curve. As a result of the smoothing- 
correction, the corrected SF{t) becomes remarkably flat 
at long timescales, and equals V^ai at timescales longer 
than ^ 20 days, which appears as an ‘ideal’ SF{t) where 
the related light curve is long enough that the edg e effects 
and aliasing are negligible ([Hughes et al.lll992n . From 
the smoothing-corrected SF{t), we see no clear asym¬ 
metry in the variations. 

We plot the /3(r) curves for both the light curves in 
Figure [5] The original /3(r) curve has a positive excess 
at long timescales, while the smoothing-corrected (3{t) 
curve appears approximately with |/3(t)| < 0.1 at the 
whole range. The (3{t) curves of the other three stitched 
sources are shown in Figure [31 all of which appear similar 
to Zw 229-15, and we find no evidence of asymmetry after 
smoothing-correction. 

4.2. Averaged Asymmetry Parameter of the Sample 

We have analyzed 19 Kepler AGNs, each of which has 
3 ^ 13 observational quarters. For each AGN, we de¬ 
rive ,0(t) curve for each of the observational quarter. 


0.01 0.10 1.00 10.00 100.00 1000.00 



Figures. /3(r) of stitched light curves of Zw 229-15. The one from 
original data is shown in grey and that after smoothing-correction 
in red. The Icr and 3 (j uncertainty ranges of 0{t), measured 
from simulated light curves assuming the variations are intrinsi¬ 
cally symmetric, i.e., with mean /3(r) = 0, are shown in filled areas 
and narrow lines with light-color, respectively. The dotted green 
line marks 0{t) = diO.l. 
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Figure 6. /3(t) of stitched light curves of W2 1925-|-50(top), W2 
1925-|-50(middle), and CGRaBS J1918-|-4937(bottom). The figure 
legends are the same as Figure [5] 

and average them to get /3(t) for each source. There 
is an alternative process, in which we average struc_ture 
funct ions of different quarters first and then derive /3(t) 
from SF{t) with Equation [2] However, with the second 
approach, the averaged asymmetry could be dominated 
by the quarter(s) with larger variation amplitudes, while 
with the first method the asymmetry parameter is sim¬ 
ply averaged over time. For the same reason, no weight 
is introducedjluring averaging. 

The mean /3 over timescales of 1 ^ 5 days and 5 ^ 20 
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Table 1 

The time-mean of all AGNs. 


Source Name 

Kepler ID 


RA 

Dec 


z 

n(g)d 

ori) 

ginal 

/^5~20 

smoothing-corrected 

type 

Zw 229-15“ 

6932990 

19 

05 

26.0 

-1-42 

27 

40 

0.028 

13 

0.04±0.07 

O.lOiO.ll 

0.01±0.05 

-0.01±0.05 

Syl 

W2 1925-1-50'' 

12158940 

19 

25 

02.2 

-1-50 

43 

14 

0.067 

11 

0.04±0.10 

-0.04±0.12 

0.02±0.06 

-0.02±0.06 

Syl 

W2R 1858-1-48 

11178007 

18 

58 

01.1 

-1-48 

50 

23 

0.079 

8 

-0.23±0.09 

-0.24±0.12 

-0.10±0.06 

-0.01±0.07 

Syl 

W2R 1904-1-37'' 

2694186 

19 

04 

58.7 

-1-37 

55 

41 

0.089 

10 

0.03±0.08 

0.09±0.13 

-0.04±0.04 

-0.04±0.06 

Syl 

W2R 1914-1-42 

6595745 

19 

14 

15.5 

-1-42 

04 

59 

0.502 

6 

-0.03±0.08 

0.06±0.15 

0.00±0.04 

0.09±0.07 

QSO‘' 

W2R 1920-1-38 

3337670 

19 

20 

47.7 

-1-38 

26 

41 

0.368 

6 

O.lOiO.lO 

0.18±0.17 

0.00±0.04 

0.00±0.06 

QSO‘' 

W2R 1931-1-43 

7610713 

19 

31 

12.5 

-1-43 

13 

27 

0.439 

7 

0.00±0.12 

O.OOiO.18 

-0.03±0.05 

-0.04±0.07 

QSO^'’" 

W2R 1910-1-38 

2837332 

19 

10 

02.5 

-1-38 

00 

09 

0.130 

6 

0.04±0.10 

0.08±0.16 

-0.05±0.06 

-0.04±0.08 

SyU 

W2R 1853-1-40 

5597763 

18 

53 

19.2 

-1-40 

53 

36 

0.625 

5 

0.01±0.12 

-0.07±0.19 

0.01±0.05 

-0.04±0.07 

QSO‘> 

W2R 1845-1-48 

10841941 

18 

45 

59.5 

-1-48 

16 

47 

0.152 

6 

0.05±0.08 

0.11±0.15 

0.00±0.05 

-0.06±0.07 

QSO‘> 

W2R 1931-1-38 

3347632 

19 

31 

15.4 

-1-38 

28 

17 

0.158 

6 

O.lliO.ll 

0.02±0.17 

0.14±0.06 

0.08±0.07 

SyU 

W2R 1926-1-42 

6690887 

19 

26 

31.0 

-1-42 

09 

59 

0.154 

6 

0.00±0.04 

-0.03±0.07 

0.00±0.03 

0.04±0.05 

BL Lac 

KA 1915-1-41 

5781475 

19 

15 

09.1 

-t41 

02 

39 

0.220 

3 

-0.27±0.20 

-0.37±0.28 

O.OOiO.lO 

-0.13±0.12 

QSO‘' 

KA 1922-1-45 

9215110 

19 

22 

11.2 

-1-45 

38 

06 

0.115 

7 

-0.07±0.06 

0.09±0.11 

-0.05±0.04 

0.01±0.06 

Syl.9 

IRXS J192949.7-I-462231 

9650715 

19 

29 

50.5 

-1-46 

22 

24 

0.127 

4 

-0.10±0.16 

-0.14±0.21 

0.09±0.09 

0.04±0.09 

Syl 

MG4 J192325-I-4754 

10663134 

19 

23 

27.2 

-1-47 

54 

17 

1.520 

11 

0.01±0.04 

0.06±0.09 

0.00±0.02 

0.02±0.04 

QSO 

MG4 J190945-I-4833 

11021406 

19 

09 

46.5 

-1-48 

34 

32 

0.513 

8 

-0.15±0.05 

-0.32±0.11 

-O.OliO.Ol 

-0.06±0.03 

QSO 

CGRaBS J1918-I-4937'' 

11606854 

19 

18 

45.6 

-1-49 

37 

55 

0.926 

11 

0.04±0.06 

0.14±0.11 

-0.02±0.03 

-0.10±0.05 

QSO 

[HB89] 1924-1-507 

12208602 

19 

26 

06.3 

-1-50 

52 

57 

1.098 

11 

-0.15±0.03 

-0.35±0.08 

O.OliO.Ol 

0.00±0.03 

Syl.5 


^ The four AGNs have stitched light curves, which have been analyzed in Section l4.ll 
^ The eight AGNs are classified as QSOs or Seyfert galaxies according to their Mj. 

^ lEdelson et al.l 1120131') classified W2R 1931+43 as a Seyfert 1 galaxy. 

^ The column shows number of observational quarters for each Kepler AGN (the same in Table (2]). 


days are ako derive^ and listed in Table [TJ The uncer¬ 
tainties of /3 i^ 5 and /? 5^20 individual sources are also 
calculated through simulations. Again, after smoothing- 
correction, most of the sources show very small asym¬ 
metry parameter (both positive and negative values are 
seen with |/3| < 0.1), consistent with zero within the sta¬ 
tistical uncertainties. For the whole sample there is not a 
general trend towards a positive or negative asymmetry. 
This indicates there is no or at most very weak variability 
asymmetry in Kepler AGNs. 

To derive better constraint on the asymmetry parame¬ 
ter, we average /^t) from all quarters of all sources. The 
uncertainties in_/3(T) are also derived from simulationfQ. 
The co-added /3 (t) averaged over all sources and its un¬ 
certainty are plotted in Figure [71 There is little differ¬ 
ence between the original and smoothing-corrected /3('r) 
curves in FigurejTl which implies that the bias due to long 
term variations with limited duration of light curves has 
been extensively reduced after averaging a large num¬ 
ber (145) of quartern. However, we note the smoothing- 
corrected co-added /3(r) has considerably smaller scatter 
(Figure [7]), as the scatter due to long term variations has 
been reduced. 

The /1 (t) averaged over timescales of 1 ~ 5 days and 
5 ~ 20 days of the whole sample and two subsamples 
(quasars and Seyfert galaxies) are listed in Table |2j We 
see that the values of the co-added /3 (t) on timescales of 
1 ~ 20 days are all consistent with zero within the small 
statistical uncertainties (again derived with simulations), 
suggesting the variations are highly symmetric. 

^ Since there is a large number (145) of Kepler quarters (real 
data), we can measure the intrinsic scatter of /3 (t) from different 
quarters and de_rive the uncertainties in the averaged /3 (t). The 
uncertainties in /3 (t) derived with this approach are slightly smaller 
than those from simulations, by a factor of 1.5 — 2.0. In this work, 
we adopt the more conservative measurements from simulations. 
Nevertheless, our conclusions are not affected by the selection. 


0.01 0.10 1.00 10.00 100.00 



0.01 0.10 1.00 10.00 100.00 
Interval (day) 

Figure 7. The average /3(r) of the whole sample including all 
sources. The figure legends are the same as Figure [5] 

5. DISGUSSION 

Using all Kepler light curves of 19 AGNs, we find the 
variability asymmetry at timescales of 1 20 days is 

rather weak. The averaged is 0.00±0.03 and -0.02±0.04 
over timescales of lr^5 days and 5~20 days, respectively, 
statistically consistent with zero. 

It’s convenient to interpret the asymmetry parame¬ 
ter in terms of shot-noise model, in which the varia¬ 
tions in AGNs are attributed to the stoch astic superpo¬ 
sition of independent discrete flares (e.g., iNegoro et al.l 
119951) . Such model provides a mathematical framework 
for a series of phy sical models (jCid Fernandes et al.ll200fll 
iFavre et al.ll2005h. inclu d ing starburst an d micro-lensing 
([Kawaguchi et al.T 119981 [Hawkins! I2002D . Generally, a 
positive asymmetry parameter /3 indicates the flare rises 
rapidly and decays gradually, such as expected in the 
starburst model, which attributes the variations in AGNs 
to random superposition of supernovae in the nuclear 
starburst region. Clearly, our detection of highly sym¬ 
metric variations in Kepler AGNs does not favor the 
starburst model. Micro-lensing m odel does predict sym¬ 
metric variations ([Hawkins! I2002D . however, for our low 
redshift sample, the probability of micro-lensing would 
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Table 2 

The time-mean I3^ of mixed samples. 


Sample 


n{q) 


original 


/ 35 ~ 


/3i, 


smoothing-corrected 
'5 ^ 5~20 


Quasars 

63 

-O.0l±0.06 

6.0O±O.l2 

-O.01±O.03 

-0.03L0.05 

Seyfert Galaxies 

76 

-0.03±0.06 

-0.04±0.10 

0.00±0.04 

-0.01±0.05 

All Sources 

145 

-0.02±0.06 

-0.02±0.10 

0.00±0.03 

-0.02±0.04 


be insignificant (jHawkinsI l2002t) . Actually, independent 
to variability asymmetry studies, these two models (star- 
burst and micro-lensing) are disfavored by observations 
which show the variations of emission lines in AGNs 
closely correlate with but lag continuum variations (e.g., 
iPeterson et alJ[2003l . 

In the disk i nstabi lity (hereafter DI) model of 
I Kawaguchi et all (|1998fl . the optical variability is as¬ 
cribed to instabilities of accretion disk as matter flows, 
and the asymmetry is due to single large scale avalanches. 
By considering a disk atmosphere emitting a power- 
law X-ray spectrum that fluctuating in time and as¬ 
suming opt ical variations simply fo llows X-ray with lit¬ 
tle delay. iKawaguchi et al.l (1199811 made Monte-Carlo 
simulations ado p ting t he cellular-automaton model of 
IMineshige et al.l (119941 ) and treated the atmosphere as 
advection-dominated accretion flow. The simulated light 
curve shows an negative asymmetry with slow rise and 
rapid decline on timesc ales of one to several h undred 
days. As simulated bv IKawaguchi et ^ (|1998ll . in Dl 
model, /3(t) = —0.1 corresponds approximately to the 
ratio of diffusion mass to inflow mass of 0.1 ^ 0.5 (cf. 
their Fig. 7), with the ratio of outer to inner disk radii 
fixed to be 20. Our strong constraints to the asymmetry 
parameter from the co-added sample therefore request an 
even higher ratio of diffusion mass to inflow mass based 
on this model. 

Theoretical calculations on the variability asymmetry 
parameter for more specific physical models are required 
to make comparison with observations. Nevertheless, the 
result of this work indicates the variations in AGNs rise 
and decay highly symmetrically. If we attribute AGN 
variations to perturbations in the accretion disk, such 
perturbations need also behave symmetrically in both 
directions of time. 

6. CONGLUSIONS 

In this paper we use the high quality light curves from 
Kepler space telescope to analyze the variability asym¬ 
metry of 19 AGNs. An asymmetry parameter /3(t) is 
introduced for quantitive description. We perform ex¬ 
tensive Monte-Carlo simulations to derive the statisti¬ 
cal uncertainties in /3(t), which can not be obtained 
through the standard error analyses approach. After 
correction of observational bias due to long term trend 
in the light curves, we find no evidence of asymmetry 
in individual sources at timescales below 20 days. For 
the whole sample there is not a general trend towards 
a positive or negative asymmetry. Co-adding data for 
all 19 AGNs, we derive an averaged /3 of 0.00±0.03 and 
-0.02±0.04 over timescales of 1^5 days and 5~20 days, 
respectively, statistically consistent with zero. Quasars 
and Seyfert galaxies tend to show similar asymmetry pa¬ 
rameter at the observed timescales. The constraint on 
longer timescale variability asymmetry is weaker as it 


requires much larger samples or much longer light curves 
which could better sample the long term variations. The 
fact that the short term optical variations in quasars and 
Seyfert galaxies are highly symmetric could put indepen¬ 
dent constraints on physical models of AGN variations. 
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